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A cluster Monte Carlo method for systems of classical spins 
with purely dipolar couplings is presented. It is tested and 
applied for finite arrays of perpendicular Ising dipoles on the 
triangular lattice. This model is a modification with long- 
range interactions of the geometrically frustrated Ising anti- 
ferromagnet. From measurements of integrated autocorrela- 
tion times for energy, magnetization and staggered magne- 
tizations, a high efficiency of the cluster Monte Carlo (MC) 
algorithm compared to a single-spin-flip algorithm is found. 
For the investigated model, a finite temperature transition is 
found which is characterized by a peak in the specific heat 
and in the staggered susceptibilities. 
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The progress in the fabrication of very small patterns 
of ferromagnetic materialsErtll opens the fascinating pos- 
sibility to study experimentally statistical properties of 
some mesoscopic or exotic magnetic systems. Dipolar 
couplingSpbetween such particles may induce magnetic 
orderingau. These couplings are often only viewed as 
additional effective anisotropy barriers for the super- 
paramagnetic relaxation of individual particles. How- 
ever, such assemblies of ferromagnetic nano-particles are 
coupled magnetic systems. They should display mag- 
netic phase-transitions for systems with mesoscopic size 
of individual magnetic moments. Recently, Wirth ct 
al. studied the properties of nanoscale ferromagnetic 
columns which were gxpwn on substrates in a scan- 
ning tunnel-microscopeo. Magnetic-force microscopic 
pictures of densely packed columns on a triangular lat- 
tice showed typical antiferromagnetic correlations where 
rows of nearest-neighbour moments point alternatingly 
up and down. 

This paper presents a numerical study on the expected 
thermodynamics for such assemblies of perpendicular 
dipoles in two dimensions. Its second aim is the intro- 
duction of a cluster Monte Carlo (MC) method for dipo- 
larly coupled systems. For the geometrically frustrated 
triangular lattice, this method allows efficient statisti- 
cal sampling down to relatively low temperatures. The 
coupled system we want to study consists of i = 1 . . . N 
Ising-spins, o~i G { — 1,1} 

1 ^-V 

H=- 2^ v{R lj )G l a j , (1) 

i,i 

where v(R) — g R~ 3 are interactions between three- 
dimensional dipoles describing the coupling between the 



perpendicular moments of the particles. All the couplings 
are antiferromagnetic. Length is measured in nearest- 
neighbour distances. If we restrict the interactions in 
(|l|) to nearest-neighbour terms we get an Ising antiferro- 
magnet. The triangular Ising antiferromagnet (TIAFM) 
is one of the few exactly solved systems with geometri- 
cal frustrations. It was shown that no long-range order 
is established at finite temperature, but the system be- 
comes critical at T = 0. However, additional interac- 
tions may lift the ground-state degeneracy and may in- 
duce phase-transitions in such underconstrained systems. 
Hence, model (Q) is interesting in its owjijjight. 

Cluster-MC methods for spin modelsETQ are only effi- 
cient in critical regions if the clusters do not percolate 
already above the critical temperature. The energy of 
model (Q) may be decomposed into two parts 

H = Y,Va + 1 X>(i2ii), (2) 

A (i,j)0NN 

where the sum of Va denotes sums of nearest-neighbour 
interactions in one of the two types of triangular plaque- 
ttes (and, with open boundary conditions, some nearest- 
neighbour interactions along the boundaries also labeled 
by Va)- This first part is exactly the TIAFM. The second 
sum contains the remaining bonds with longer range. 

Clusters for the TIAFM can be defined by sampling 
the sums of energies corresponding to elementary trian- 
gular plaquetteso. These clusters are the correct critical 
clusters of the TIAFM percolating at T = 0. Applying 
this cluster definition to the first part of ^ the dipolar 
system is mapped on a system of clusters which interact 
through the long-range interactions generated from terms 
in the second sum of (^). For flipping clusters the usual 
Metropolis algorithm is used. The algorithm of alternat- 
ingly mapping (1) onto interacting clusters via (2) and 
sampling this mapped system is equivalent to sampling 
the original system. 

We have to show ergodicity and detailed balance to jus- 
tify the method. The procedure of cluster decomposition 
may result in deleting all bonds which maps the origi- 
nal system onto itself. Thus, the occurrence of single- 
spin flips ensuring ergodicity can be checked a posteriori. 
Detailed balance may be demonstrated by the prescrip- 
tion for general cluster algorithms given by Kandel and 
DomanyQ which shows that bonds may be deleted, frozen, 
or left unchanged to generate systems of interacting clus- 
ters equivalent to the original system. 

The method is applied for finite arrays of hexago- 
nal shape with N = 61 . . . 2977 spins. The follow- 
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ing measure for antiferromagnetic correlations is used: 
Spins are labeled cr M „ for sites given by + Wi with 
nearest-neighbour vectors = (±v / 3/2, 1/2). Three 
types of staggered magnetizations may be calculated by 

™ J = i/JV ^(-i)", m 11 = i/N ^(-i)", 

and m 111 = 1/TV (^(-l)^". For the hexagonal 
arrays, the three types are equivalent. Thus we use the 
average m a = 1/3 (ra 1 + m + to 111 ). 

Here, our algorithm is first used in the Swendsen- 
Wang (SW) form by deleting/freezing bonds for all pla- 
quettes of a randomly chosen type. It turns out that 
after mapping onto the system of interacting clusters 
best performance was achieved by attempts to flip just 
one (randomly chosen or the largest) cluster. Integrated 
auto-correlation times for energy, magnetization, and 
for the staggered magnetizations to 1 - 111 were calculated 
as measure of efficiency by a self-consistent windowing 
methocO. Independent isothermal runs were started 
from slowly cooled and relaxed configurations. For sizes 
TV = 61 ... 469 between several 10 4 and 10 5 cluster-MC 
sweeps (MCS, i.e., TV attempted cluster moves), for the 
two larger sizes only 10 3 ' 4 MCS could be performed. 
As expected, the autocorrelations for the magnetization 
decay rapidly at all temperatures. The autocorrelations 
for energy are similar to those of the staggered magneti- 
zations. Thus, we only show the average integrated au- 
tocorrelation time r a for the staggered magnetizations in 
Fig. 1(a) and compare it to runs with the single-spin- flip 
Metropolis algorithm for the two smallest system sizes. 
Fig. 1(b) shows the average size <c> of flipped clusters. 
The strong increase of < c > towards low temperature 
means that all spins are increasingly correlated. This 
transition becomes sharper with larger TV. The size dis- 
tribution of flipped clusters always consists of a rapidly 
decreasing part starting with single-spin flips. For low 
temperatures and small system sizes clusters with size 
close to TV are found also. 

Still, the corresponding r a does not grow dramati- 
cally. For TV > 469 we could measure r Q only down 
to T/g ~ 0.19, where it becomes essentially indepen- 
dent of TV while <c>~ 2. The single-spin- flip algorithm 
yields far larger r a which in a critical region should grow 
with system size as r a cx TV 2 / 2 with z ~ 2. This com- 
parison demonstrates the efficiency of this cluster-MC. 
However, the numerical cost of the SW-algorithm scales 
as TVr a . Therefore, a Wolff- type algorithm!] was imple- 
mented for further production runs shown below. Only a 
single cluster is grown starting from a randomly chosen 
plaquette. The performance of this algorithm is similar 
to the SW-algorithm but the numerical effort is smaller. 
Its numerical cost scales as < c > r a . 

Cluster-MC runs with several 10 5 MCS were performed 
for TV < 721. For TV > 1519 serious problems in the 
critical region T/g < 0.18 arise. Cycling temperature 
displays hysteresis of the energy. Measurements could 
be done for runs with only few 10 3 MCS in the critical 
region with TV > 1519. Here, configurations both from 



cooling and heating were first relaxed for longer times be- 
fore start of measurements to ensure that they are close to 
equilibrium. Sampling was performed for energy per spin 
E, specific heat C = N{<E 2 > - <E> 2 )/T 2 , the stag- 
gered magnetizations, the corresponding susceptibilities 
X x = TV(< (m x ) 2 > - < m x > 2 )/T and the fourth-order 
Binder cumulant U x = 1- < (m x ) 4 > /(3 < (m x ) 2 > 2 ) 
(X = I, II, III). Using the histogram method by Ferren- 
berg and Swendsenliil these various quantities were eval- 
uated as functions of temperature in the critical region. 
Results are shown in Figs. 2 and 3. The diverging peak 
in the specific heat (Fig. 2(b)) indicates the presence 
of a transition at or closely below T/g ~ 0.18. Simi- 
larly, the peak for the average staggered susceptibility 
Xa shows that the transition (Fig. 3(a)) is related to an- 
tiferromagnetic ordering. Near T/g ~ 0.18, m a steeply 
increases reaching values ~ 0.3 to 0.4 (not shown here). 
However, the average Binder-cumulant U a for the various 
system sizes does not provide a clear indication for a con- 
tinuous transition by a unique crossing of Z7 Q (T)-curves 
for the different system sizes at one critical temperature 
(Fig. 3(b)). The transition might be first order instead as 
indicated by the observed hysteresis when cycling tem- 
perature. Then the apparent correlation of system sizes 
TV < 1000 corresponds only to the size of typical nuclei 
of the low-temperature phase. On the other hand, the 
behaviour of the U a (T)-curves which seem to merge be- 
low T/g — 0.18, is consistent with usual observations at 
a Berezinskii-Kosterlitz-Thouless-transition and a line of 
critical points at low temperature. Possibly, the peak in 
the specific heat is only a precursor of this transition, or 
two transitions are present. 

In conclusion a cluster-MC method is presented which 
may be used for dipolarly coupled spin-systems. It works 
by generating interacting clusters from sampling of short- 
range interactions. It is shown to be efficient for the 
specific case of Ising dipoles on the triangular lattice. 
This success relies on the ability to define clusters for this 
geometrically frustrated system which percolate only at 
T = 0. However, the method can be generalized to vari- 
ous other systems with long-range interactions, including 
systems with continuous spinal The investigation of the 
exact nature of the finite-temperature transition found 
with this MC-method remains a task for the future. 
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FIG. 1. (a) Average integrated autocorrelation time for 
staggered magnetizations (see text) versus temperature. Sym- 
bols correspond to independent isothermal runs. Comparison 
between Metropolis single-spin-flip and cluster algorithm. For 
the largest size N = 2977, only one run at T/g = 0.190 is 
shown (big circle), (b) Corresponding average cluster sizes 
< c > for different system sizes. 
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FIG. 2. (a) Energy per spin (b) specific heat in the region 
of the transition temperature of Ising dipoles on hexagonal 
triangular-lattice arrays with various sizes N. 
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FIG. 3. (a) Average staggered susceptibility \a = 1/3 (x^ 
+X^)- (b) Corresponding average fourth-order Binder 
cumulant in the critical temperature region. 
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